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Factors  associated  with  short  circuit  current  density  (Jsc)  and  open  circuit  photovoltage  (\70C)  of  dye 
sensitized  solar  cells  (DSSCs)  have  been  analyzed  through  DFT  and  TDDFT  calculations  to  explore  the 
origin  of  the  significant  performance  differences  with  only  tiny  structure  difference  (1.24%  for  1  and 
8.21%  for  2)  (Advanced  Functional  Materials  2012,  22, 1291—1302).  Our  results  reveal  that  the  insertion  of 
phenyl  ring  in  2  enlarges  the  distance  between  the  dye  cation  hole  and  the  semiconductor  surface  and 
makes  the  benzothiadiazole  (BTDA)  unit,  which  has  strong  interaction  with  the  electrolyte,  far  away  from 
the  semiconductor,  resulting  in  a  decreased  charge  recombination  rate  compared  with  that  of  1.  How¬ 
ever,  the  insertion  of  phenyl  ring  also  results  in  a  distortion  of  the  molecular  structure,  leading  to  a 
decreased  light  harvesting  ability.  Hence,  two  dyes  (6  and  7)  derived  from  2  with  better  conjugation 
degree,  farther  position  of  BTDA  unit  and  longer  molecular  length  have  been  designed  to  keep  the  ad¬ 
vantages  and  overcome  the  disadvantages  of  2  simultaneously.  The  results  demonstrate  that  we  get  the 
desired  properties  of  dyes  through  reasonable  molecular  design,  and  these  two  dyes  could  be  promising 
candidates  in  DSSC  field  and  further  improve  the  performance  of  the  cell. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 
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As  one  of  the  most  promising  renewable  energy  devices,  dye 
sensitized  solar  cells  (DSSCs)  are  receiving  tremendous  attention 
due  to  their  low  production  costs,  good  flexibility,  light  weight  and 
easy  production  processes  [1—6].  Up  to  now,  the  overall  power 
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conversion  efficiency  (PCE,  77)  of  DSSCs  has  reached  13.0%  [7  re¬ 
ported  by  Gratzel  and  co-workers  [8],  however,  which  is  still  lower 
than  that  of  the  traditional  Si-based  solar  cells  [9  .  Thus,  lots  of 
researches  are  still  devoted  to  improving  the  PCE  of  DSSCs 
[1,8,10-16].  Generally,  a  typical  DSSC  is  mainly  composed  of 
semiconductor  material,  dye  sensitizer  and  electrolyte.  There  is  no 
doubt  that  the  development  of  new  efficient  sensitizers  is  an 
effective  way  to  enhance  the  efficiency  of  DSSCs  in  view  of  the  key 
roles  of  the  sensitizer  played  in  the  cell,  such  as  affecting  the 
lighting  harvesting  efficiency,  the  charge  transfer  process  and  ul¬ 
timately  the  PCE  of  the  cell. 

Currently,  there  are  mainly  three  kind  of  dyes  commonly  used  in 
DSSC  [2],  including  ruthenium  dyes  [16-18  ,  porphyrin  dyes  [19,20] 
and  metal-free  organic  dyes  [21-23  .  Compared  with  the  other  two 
types  of  dyes,  organic  D-tu-A  dyes  have  received  an  ever-increasing 
attention  due  to  their  lower  cost,  easy  of  fabrication,  relative  higher 
molar  extinction  coefficients  and  environmental  friendly.  Impres¬ 
sively,  the  efficiency  of  DSSC  based  on  champion  organic  dye  (C219) 
is  up  to  10.3%  [24  .  However,  organic  dyes  also  have  their  own 
drawbacks  such  as  the  narrower  absorption  spectrum,  and  the  faster 


charge  recombination  due  to  the  presence  of  the  heteroatom  [25]. 
Thus,  great  efforts  have  been  devoted  to  the  design  and  optimization 
of  novel  dyes  to  overcome  these  shortcomings  [26-33  .  For 
example,  Haid  et  al.  synthesized  two  triphenylamine-based  (TPA) 
organic  dyes  [34  with  tiny  structure  difference,  as  shown  in  Fig.  1.  It 
is  found  that  the  insertion  of  phenyl  ring  between  the  benzothia- 
diazole  (BTDA)  unit  and  the  cyanoacrylic  acid  group  of  2  can  lead  to  a 
6.5  fold  increase  in  cell  efficiency  (1.24%  for  1  and  8.21%  for  2)  [34]. 
How  could  this  happen?  In  this  paper,  our  initial  interest  focused  on 
exploring  the  physical  origin  of  the  performance  difference  based  on 
two  organic  dyes  only  with  tiny  structure  difference.  Here,  in  order 
to  shed  light  on  the  reason  of  the  performance  difference,  key  pa¬ 
rameters  including  open  circuit  photovoltage  (Uoc)  and  short  circuit 
current  density  (Jsc )  were  systematically  investigated  in  detail  based 
on  density  functional  theory  (DFT)  and  time-dependent  DFT 
(TDDFT)  calculations.  In  addition,  we  also  designed  two  dyes  (6  and 
7)  with  stronger  light  harvesting  ability  and  slower  charge  recom¬ 
bination,  which  will  be  potential  sensitizers  for  DSSCs.  We  hope  our 
work  can  provide  guidance  to  the  design  of  new  dyes  and  will  move 
the  DSSCs  field  forward. 
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Fig.  1.  The  structures  of  studied  organic  dyes. 
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2.  Computational  details 

All  the  calculations  are  performed  in  Gaussian09  code  [35  at 
DFT  and  TDDFT  level.  The  ground  state  molecular  structures  of  dyes 
under  study  were  optimized  at  B3LYP/6-31G(d)  level,  followed  by 
frequency  calculations  to  make  sure  we  got  stable  structures  [36]. 
The  absorption  spectra  were  calculated  at  CPCM/CAM-B3LYP/6- 
31G(d)  level  considering  real  solvent  effect  (dichloromethane) 
with  conductor-like  polarizable  continuum  model  [37]  and  the 
long-range  correction  [38  .  To  evaluate  the  extent  of  charge 
recombination  process  between  dyes  and  electrolyte,  we  also 
calculated  the  interaction  energy  between  dye  and  I2  with  M06-2X 
functional  [39-41  ,  which  is  considered  as  an  suitable  functional 
for  the  description  of  weak  interaction,  together  with  6-31  G(d) 
basis  set  for  C,  H,  N,  O,  S  and  LANL2DZ  basis  set  and  corresponding 
effective  core  potential  for  I  atom.  And,  the  counterpoise  method 
was  used  to  compensate  the  basis  set  superposition  error  (BSSE) 
[42,43]. 

For  the  semiconductor  surface  effects,  we  take  (Ti02)3g  cluster 
[44]  obtained  by  cutting  an  anatase  slab  primarily  exposing  the  101 
surface  into  consideration  as  shown  in  Fig.  SI,  because  this  surface 
has  been  demonstrated  to  be  the  most  stable  and  typically 
observed  surface  for  this  material  [45,46  .  De  Angelis  and  co¬ 
workers  also  found  that  the  excitation  energy  of  (Ti02)3g  cluster 
is  in  good  agreement  with  the  experimental  results  of  the  Ti02 
semiconductor  47,48].  Thus  we  believe  that  it  is  reasonable  to  use 
this  model  to  simulate  the  effects  of  the  semiconductor  surface.  The 
structures  of  dyes  binding  to  the  (Ti02)3g  cluster  in  vacuum  was 
optimized  by  SIESTA  code  [49-51  ,  and  the  computational  details 
are  listed  in  SI  (Supporting  information).  Besides,  the  Density  of 
States  (DOS)  profiles  used  to  estimate  the  change  of  conduction 
band  energy  level  was  obtained  at  C-PCM/B3LYP/SVP  level  [52  . 

3.  Results  and  discussion 

The  PCE  of  the  cell  has  a  close  relationship  with  the  Jsc,  Voc  and 
the  fill  factor  (FF).  And  it  can  be  estimated  as  [53]: 


the  unit  charge.  It  is  obvious  that  ACB  and  nc  are  the  factors  which 
could  change  the  Voc  of  the  cell.  In  fact,  a  large  number  of  theo¬ 
retical  and  experimental  investigations  [55-60]  have  demon¬ 
strated  that  for  the  cells  under  the  same  experimental  conditions 
with  just  different  sensitizers,  the  main  reason  for  the  Voc  differ¬ 
ence  is  the  different  ACB  and  different  extent  of  charge  recombi¬ 
nation  affecting  nc. 

Actually  there  are  many  factors  could  induce  ACB,  such  as,  sol¬ 
vent,  ions  in  the  electrolyte  and  dye  sensitizers.  The  express  is 
presented  as  follows  [56,58,59]: 

ACBtot  =  ACBsoiv  +  ACBjons  +  ACB^ye  (3) 


ACBtot  means  total  Ti02  conduction  band  shift,  ACBsoiv,  ACBions  and 
ACBdye  represent  the  shifts  induced  by  the  solvent,  ions  and 
adsorbed  sensitizer  in  electrolyte,  respectively.  However,  under  the 
same  experimental  condition,  ACBsoiv  and  ACBions  could  be 
considered  as  the  constant.  Thus  it  is  easy  to  get  a  conclusion  that 
the  conduction  band  energy  shift  of  dye  sensitizer  has  a  significant 
effect  on  the  ACBtot  and  finally  have  an  important  influence  on  Voc 
based  on  Eq.  (2). 

For  another  important  factor  affecting  the  performance  of  the 
cell  from  the  sensitizers  point  of  view,  it  is  the  charge  recombina¬ 
tion  as  described  above,  which  could  affect  nc.  There  are  mainly  two 
charge  recombination  processes  in  the  solar  cell,  one  is  the 
recombination  between  the  injected  electrons  and  dye  cation 
(fcreci),  and  the  other  is  that  between  the  injected  electrons  and 
electrolyte  (krec 2)-  Detailed  descriptions  are  listed  as  follows: 

/<rec  =  fcrecl  +  ^rec2  (4) 

Based  on  the  electron-transfer  theory  61-63],  /creci  could  be 
expressed  as  64] : 
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where  Pmax  and  P[n  are  the  maximum  and  incident  power  of  the 
DSSCs,  respectively. 

As  discussed  in  the  part  of  the  introduction,  we  know  that  1  and 
2  synthesized  by  Haid  et  al.  34]  with  tiny  structure  difference  have 
significant  different  photovoltaic  performance  data  C/Sci  - 
3.40  mA  cm-2,  Voci:  489  mV,  FFi:  0.74;  Js c2:  18.47  mA  cm“2,  V0C2- 
640  mV,  FF2:  0.69).  Thus  in  order  to  get  a  deep  understand  on  the 
essential  reasons  of  the  significant  difference  in  77  of  the  solar  cells 
based  on  these  two  dyes  from  the  theoretical  point  of  view,  we  will 
analyze  how  the  molecular  structure  of  the  sensitizer  affects  the  Jsc 
and  Voc. 


3.1.  Theoretical  explanations  for  the  Voc  difference 
For  Voc  in  DSSCs,  it  can  be  determined  by  [54]: 

yoc  =  £cb  +  ACB  +  — lnf'Y)  -  (2) 

q  q  \NC )  q 

here,  Fcb  is  the  conduction  band  energy  level,  ACB  is  the  shift  of  the 
conduction  band  energy  level.  nc  is  the  number  of  the  electrons  in 
the  conduction  band  of  semiconductor,  Nc  is  the  electron  density  of 
state,  Fredox  is  the  oxidation  potential  of  the  redox  couple  in  the 
electrolyte,  kg  is  the  Boltzmann  constant,  T  is  the  temperature,  q  is 


b,  AG°  are  the  electron  coupling  and  the  free  energy  between 
the  donor  and  acceptor  states,  respectively,  and  A  is  the  reorgani¬ 
zation  energy.  For  theH^B,  it  is  related  to  electron  tunnelling 
through  a  potential  barrier  (/?,  which  is  a  function  of  the  barrier 
height  of  intervening  media)  and  shows  an  exponential  increase 
with  the  increasing  spatial  distance  (r)  between  donor  and  acceptor 
states.  Obviously,  /<reci  is  determined  by  AG°,  A,  {3  and  r.  Among  the 
four  parameters,  A  is  significantly  controlled  by  the  solvent,  and  (3  is 
associated  with  the  intervening  media,  which  means  that  they  are 
not  molecular  control.  Assuming  all  the  DSSCs  device  are  fabricated 
in  same  experiment  condition,  the  A  and  (3  would  be  same.  In 
addition,  for  AG°,  which  is  defined  as  the  difference  between  the 
ground  state  oxidation  potential  of  dyes  and  the  Fermi  level  of 
semiconductor,  Clifford  et  al.  has  demonstrated  that  /<reci  is  not 
sensitive  to  the  AG°  [64  .  Thus,  we  focus  on  the  spatial  distance  r 
between  donor  and  acceptor  to  estimate  /creci-  Here,  donor  is 
defined  as  the  injected  electron  in  the  conduction  band  of  semi¬ 
conductor  and  acceptor  as  the  dye  cation  hole  (dye  cation  HOMO  (3 
orbital  was  used). 

For  the  kreC2,  as  we  all  know,  the  higher  concentration  of  reac¬ 
tant  is,  the  higher  reaction  rate  is.  Some  experiments  have  revealed 
that  some  atoms  or  groups  could  provide  binding  sites  for  iodine 
through  halogen  bond  [65-67],  which  would  increase  the  local 
concentration  of  iodine  in  the  vicinity  of  the  semiconductor  surface 
[60,68-72].  And  De  Angelis  and  co-workers  have  demonstrated 
that  the  more  predominant  interaction  sites  the  dye  has  and  the 
shorter  distance  between  iodine  and  semiconductor  surface  is,  the 
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larger  the  /<rec 2  is  [25,60].  And  this  strategy  has  been  widely  applied 
in  theoretical  research  of  DSSCs  [57,58,73  .  Thus,  we  calculated  the 
interaction  energy  between  dyes  and  iodine  and  the  approximate 
distance  between  iodine  and  semiconductor  surface  to  analyze  /<rec 2 
qualitatively. 

3  A  A.  Molecular  structure  effects  on  the  conduction  band  energy 
level  shift,  ACB 

In  order  to  get  the  conduction  band  energy  shift  induced  by 
the  adsorption  of  sensitizer,  we  firstly  carried  out  full  geometry 
optimization  of  dyes  adsorbed  onto  Ti02  surface  using  the 
preferred  bidentate  bridging  mode  with  deprotonation  after 
adsorption  [74,75  .  The  optimized  structures  of  dye/(Ti02)38  are 
shown  in  Fig.  2.  Then,  based  on  the  optimized  geometries,  the 
DOS  profiles  for  Ti02  with  dyes  and  Partial  Density  of  States 
(PDOS)  profiles  for  pure  Ti02  at  the  B3LYP/SVP  level  in  acetonitrile 
solvent  were  analyzed  [56].  As  shown  in  Fig.  3,  the  corresponding 
ACB  for  1  and  2  are  0.263  eV  and  0.256  eV,  suggesting  that  1  and 
2  have  similar  shift  in  conduction  band  edge.  What  is  the  origin  of 
conduction  band  shift?  Actually,  there  are  two  effects  proposed 
by  Ronca  et  al.  [56]  when  dye  adsorbs  onto  the  semiconductor. 
One  is  the  charge  transfer  effect  due  to  electron  density 


1-Ti02 


arrangement,  and  the  other  is  the  electrostatic  effect  due  to  the 
dipole  moment.  Thus,  we  also  quantitatively  estimate  this  two 
effects  through  the  entity  of  charge  transfer  (qCT)  originated  from 
the  charge  displacement  curve  and  the  average  electrostatic  po¬ 
tential  (VEl)  generated  by  the  molecular  charges  evaluated  on  the 
first  layer  of  the  Ti02  cluster  as  shown  in  fable  1.  The  computa¬ 
tional  details  should  be  found  in  the  paper  of  Ronca  et  al.  56]  and 
our  previous  works  [58,73  .  From  Table  1,  it  is  found  that  the  qCT 
of  1  and  2  are  0.38  and  0.36  respectively,  and  VEi  of  them  of 
all  -0.21.  Obviously,  the  similar  ACB  should  be  owing  to  the 
similar  qCT  and  VEE.  In  any  case,  it  should  be  found  that  ACB  is  not 
the  dominant  factor  resulting  in  the  significant  difference  of  Voc 
between  1  and  2.  Thus,  we  discussed  the  charge  recombination  in 
the  following  in  an  attempt  to  rationalize  the  significant  differ¬ 
ence  of  i/oc. 

3 A. 2.  Molecular  structure  effects  on  the  charge  recombination 
process,  nc 

As  mentioned  above,  nc,  which  is  closely  associated  with  charge 
recombination  including  /creci  and  krec2,  also  plays  an  important  role 
in  affecting  i/oc.  Therefore,  in  the  following  we  will  discuss  how  the 
molecular  structure  difference  affect  /<reci  and  /<rec 2- 


Fig.  2.  The  optimized  molecular  structures  of  l-(Ti02)38  and  2-(Ti02)38  by  SIESTA  code. 
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E  (eV) 

Fig.  3.  Total  and  partial  density  of  states  (DOS)  of  l-(Ti02)38  and  2-(Ti02)38. 


As  for  /<reci,  it  is  mainly  determined  by  the  hole-surface  distance 
r  between  donor  and  acceptor  states  according  to  the  nonadiabatic 
electron  transfer  theory.  Moreover,  Durrant  and  his  co-workers  [64] 
also  found  that,  in  general,  DSSCs  employing  smaller  dyes  exhibit 
faster  /creci  because  of  the  shorter  r  between  the  dye  cation  hole  and 
semiconductor  surface.  Thus,  we  followed  their  strategy  and 
calculated  the  spatial  distance  r  by  multiplying  the  percentage 
contribution  of  each  atom  to  the  HOMO  and  the  corresponding 
distance  between  the  carbon  atom  of  the  carboxylic  group  and  each 
atom  of  the  dye  cation  to  evaluate  /creci*  The  results  are  shown  in 
Fig.  4.  From  Eq.  (6),  we  know  that  the  H^B  exhibits  an  exponential 
increase  with  respect  to  the  distance  r,  the  longer  r,  the  larger  Hj 3. 
Hence,  the  3.6  A  lower  of  1  on  r  than  that  of  2  could  lead  to  sig¬ 
nificant  decrease  of  /<reci. 

To  analyze  the  extent  of  /<rec 2  of  different  dyes,  the  interaction 
energies  between  dye  and  I2  and  their  corresponding  approximated 
distances  were  calculated  and  shown  in  Fig.  5.  We  can  see  that,  for 
1,  it  has  three  predominant  sites,  1-N2  (-5.74  kcal  mol”1;  1.3  A),  1- 
N3  (-5.67  kcal  mol”1;  7.1  A)  and  1-N1  (-5.65  kcal  mol'1;  10.1  A), 
respectively,  implying  that  I2  prefer  to  surround  the  BTDA  unit  and 
CN  group  rather  than  the  others.  With  respect  to  2,  the  predomi¬ 
nant  sites  are  at  2-N2  (-6.73  kcal  mol”1;  6.4  A)  and  2-N1 
(-6.13  kcal  mol”1;  14.1  A),  also  showing  that  BTDA  unit  has  strong 
interaction  with  I2.  Compared  2  with  1,  following  the  conclusions  of 
De  Angelis  and  his  co-workers,  we  believe  that  2  has  a  smaller  /<rec 2 
than  that  of  1  because  2  has  less  predominant  sites  and  each  of  the 
sites  are  farther  than  that  of  1.  The  calculated  results  suggest  that 
BTDA  unit  could  interact  with  I2  strongly,  and  that  the  insertion  of 


phenyl  ring  could  lead  to  the  BTDA  unit  much  farther  away  from  the 
Ti02  surface,  which  results  in  the  /<rec 2  of  2  smaller  than  1.  Obvi¬ 
ously,  BTDA  units  play  important  roles  in  controlling  /<rec 2- 

It  is  worth  mentioning  that  there  are  different  kinds  of  dye- 
iodine  complexes,  such  as  dye-I2  and  I3,  and  the  dye-I2  is  the 
most  important  complex  as  its  charge  recombination  rate  is  2  order 
of  magnitude  faster  than  the  dye-Ff  complex  which  was  demon¬ 
strated  by  Green  et  al.  [76  ,  and  this  complex  could  accelerate  /<rec 2 
has  also  been  observed  in  some  related  experiments  [69,70,77  .  In 
addition,  there  are  some  other  factors  could  also  affect  the  charge 
recombination  process,  such  as  the  blocking  effects  induced  by  the 
alky  chain  on  the  sensitizer,  the  addictive  in  the  electrolyte,  the 
thickness  of  the  semiconductor  and  so  on  2,78].  In  our  manuscript, 
these  factors  were  not  taken  into  consideration.  Because  for  the 
DSSCs  with  just  different  sensitizers,  we  assume  that  they  have  the 
same  experimental  conditions,  and  we  believe  that  the  in¬ 
vestigations  for  the  inherent  charge  recombination  of  the  sensi¬ 
tizers  could  be  more  important  for  the  further  optimization  of  the 
sensitizers.  From  the  above  discussion,  taking  both  /<reci  and  /<rec 2 
into  consideration,  we  found  that  it  is  the  fast  /<rec  that  lead  to  the 
lower  Voc  °f  1  than  that  of  2,  which  is  in  agreement  with  the 
experiment  [34  .  And  in  the  next,  we  will  discuss  how  the  molec¬ 
ular  structures  induce  the  Jsc  difference. 

3.2.  Theoretical  explanations  for  the  Jsc  difference 

Jsc ,  as  another  key  parameter  determining  the  PCE  of  the  cell, 
could  be  expressed  as  [53]; 

Jsc  =  J  eIPCE(A)/s(A)dA  (7) 

I  PCE  (A)  =  LHE(A)<f’inj?7C0[[  (8) 

where  IPCE(A)  is  incident  photon-to-current  conversion  efficiency 
at  a  fixed  wavelength,  which  is  the  product  of  LHE(A)  (light  har¬ 
vesting  efficiency  at  the  certain  wavelength),  $inj  (electron  injec¬ 
tion  efficiency)  and  rjcoii  (charge  collection  efficiency).  And  /S(A)  is 
the  corresponding  photon  flux  based  on  AM  1.5  G  (global)  solar 
radiation  spectrum. 

In  terms  of  LHE(A),  it  can  be  described  by  the  following  equation 

[79,80]: 

LHE(A)  =  1  -  Kr£Wbc  (9) 

here,  e(A)  is  the  molar  absorption  coefficient  at  certain  wavelength, 
b  is  the  thickness  of  film,  and  c  is  the  concentration  of  dye.  For 
discussing  LHE(A)  conveniently,  b  and  c  are  taken  as  10  pm  and 
300  mmol  L_1,  respectively  [81,82]. 

For  $inj,  it  is  related  to  the  injection  driving  force  (AGmj)  defined 
as  the  difference  between  the  excited  state  oxidation  potential  of 
dye  and  the  conduction  band  energy  level  of  semiconductor.  Thus, 
$inj  can  be  estimated  qualitatively  as  follow  [57,83-85]: 


AGjnj  «  fcjnj  «  ^inj 

(10) 

AGinj  =  Edye*  “  £CB 

(11) 

^dye’  =  ^dye  —  ^max 

(12) 

here,  Edye*  and  Edye  are  the  excited  and  ground  states  oxidation 
potential,  respectively.  Amax  is  the  maximum  absorption  wave¬ 
length,  ECb  is  the  conduction  band  energy  level  of  semiconductor, 
and  we  adopted  the  experimental  value  -4.0  eV  as  our  standard  Eqb 
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Table  1 

Calculated  conduction  band  shift  (ACB,  eV),  entity  of  charge  transfer  ( qCJ ,  e),  and 
average  electrostatic  potential  on  the  first  layer  of  Ti02  cluster  (VEL,  eV)  for  1  and  2. 


Dye 

ACB 

CT 

Q 

^EL 

1 

0.263 

0.38 

-0.21 

2 

0.256 

0.36 

-0.21 

in  this  paper  [86].  Based  the  Eq.  (10),  we  can  estimate  the  &mj  by 
comparing  the  AGinj. 

As  for  r/c on,  it  is  not  only  related  to  the  electrode  materials,  but 
also  have  close  relationship  with  the  charge  recombination  process. 
In  an  ideal  DSSC  device,  ^con  should  be  one.  However,  under  the 
practical  operation,  ^con  can  not  reach  100%  because  of  the  charge 
recombination  process.  For  the  same  semiconductor,  the  larger  /<reo 
the  lower  ^con. 

In  addition,  it  is  worth  to  note  that  light  harvesting  is  not  only 
associated  with  LHE(A),  but  also  related  to  the  absorption  range. 
Thus,  7S(A)  need  to  be  taken  into  account  when  light  harvesting  is 
estimated  according  to  Eq.  (7)  [53,79].  If  we  assume  the  unity  ef¬ 
ficiency  for  $inj  and  77con,  the  calculated  Jsc  could  be  the  theoretical 
maximum  limit  J™ax  which  could  quantitatively  reflect  the  overall 
light  harvesting  ability  of  dyes,  in  other  words,  the  matching  degree 
between  the  absorption  spectrum  of  the  sensitizer  and  the  photon 
flux  spectrum. 

As  a  result,  it  is  obvious  that  the  overall  light  harvesting  ability 
(theoretical  maximum  limit J™ax),  $inj  and  77con  are  the  key  factors 
associated  with  the  sensitizer  affecting  Jsc,  And  thus  in  the 
following,  we  will  give  a  detailed  interpretation  for  the  difference  of 
Jsc  from  these  aspects. 

3.2.1.  Molecular  structure  effects  on  the  light  harvesting 

In  order  to  evaluate  the  overall  light  harvesting  ability  (theo¬ 
retical  maximum  limit  J™ax),  we  firstly  calculated  the  maximum 
absorption  wavelengths,  vertical  excitation  energies,  and  nature  of 
the  transitions  of  the  dyes  under  study  and  the  results  were  sum¬ 
marized  in  able  2.  From  Table  2,  we  can  see  that  the  maximum 
absorption  peak  for  1  and  2  is  at  561  and  501  nm,  respectively,  both 
corresponding  to  the  transition  of  So  to  Si,  which  can  be  assigned  as 
the  promotion  from  HOMO  to  LUMO  and  HOMO-1  to  LUMO.  And 
the  maximum  absorption  peaks  are  in  good  agreement  with  the 
experimental  values  (570  nm  and  515  nm  for  1  and  2,  respectively) 
[34  ,  validating  the  rationality  of  our  method  and  basis  set  we 
choosed.  We  also  find  that  the  insertion  of  phenyl  ring  could  cause  a 


1  2 


Fig.  4.  The  distance  between  dye  cation  hole  and  Ti02  surface  of  1  and  2. 


blue-shifted  spectrum,  resulting  from  the  distortion  of  molecular 
structure  of  2  as  shown  in  Fig.  S2.  Inspection  of  the  distribution  of 
frontier  molecular  orbitals  of  1  and  2  depicted  in  Fig.  S3  shows  that 
the  HOMO-1  s  and  HOMOs  of  both  dyes  are  mainly  located  at  TPA 
donors  and  thiophene  rings,  and  the  LUMOs  are  mainly  located  at 
BTDA  units  and  cyanoacrylic  acid  groups,  implying  that  the  excited 
state  transition  (So-Si)  is  a  charge-transfer  transition. 

To  give  an  intuitive  impression  on  light  absorption,  we  simu¬ 
lated  the  UV/Vis  absorption  spectra  of  1  and  2  shown  in  Fig.  6 
through  the  Multiwfn  code  [87,88]  at  CPCM/TD-CAM-B3LYP/6- 
31G(d)  level.  As  discussed  above,  we  know  that  the  theoretical 
maximum  limit  J™axis  a  really  important  factor  affecting  the  Jsc  of 
the  cell,  and  it  takes  the  absorption  range,  absorptance  and  the 
photon  flux  spectrum  into  consideration  at  the  same  time.  Thus, 
according  to  Eqs.  (7)-(9),  we  get  the  LHE  curves  and  their  corre¬ 
sponding  J™ax  as  shown  in  Fig.  6.  The  J™ax  of  1  and  2  are 
29.83  mA  cm-2  and  22.05  mA  cm-2,  respectively,  implying  that  1 
should  have  better  overall  light  harvesting  ability  than  that  of  2 
leading  to  a  higher  Js c.  Moreover,  from  the  LHE  curves,  it  could 
found  that  when  the  molar  absorption  coefficient  at  certain 
wavelength  reaches  some  value,  its  corresponding  LHE(A)  will  be 
close  to  one.  And  the  larger  J™ax  of  1  than  that  of  2  is  mainly  due  to 
the  red  shift  of  the  absorption.  However,  based  on  the  experimental 
value,  we  know  that  the  Jsc  of  dye  2  is  far  greater  than  that  of  dye  1. 
The  contrary  tendency  indicates  that  there  must  be  some  electronic 
loss  processes  in  the  whole  circulation  system  of  1-based  DSSCs, 
causing  that  its  Jsc  was  reduced  rapidly.  Thus,  we  will  focus  on  the 
$inj  and  rjcoi,  which  are  related  to  the  electronic  loss  in  the 
following. 

3.2.2.  Molecular  structure  effects  on  the  electron  injection  efficiency 
(@inj)  and  electron  collection  efficiency  (ricou) 

The  driving  force  was  used  to  evaluate  the  $inj,  and  Islam  et  al. 
has  demonstrated  that  when  the  injection  driving  force  is  above 
0.2  eV,  the  $inj  will  be  unit  [89  .  Table  3  shows  obviously  that  these 
two  dyes  have  enough  driving  force  to  insure  the  photoexcited 
electron  injection,  which  is  in  good  agreement  with  the  experiment 
data  (0inji:  94.8%,  0^2:  94.7%)  [34].  Thus  it  is  obvious  that  $inj  can 
not  determine  the  difference  in  Jsc  between  two  dyes.  How  does 
77Coii  impact  it? 

The  77Coii  is  mainly  determined  by  /<reo  leading  to  the  loosing  of 
injected  electrons.  From  the  Section  3.1.2,  we  know  that  the  hole- 
surface  distance  r  of  2  is  larger  and  its  predominant  interaction 
site  with  I2  is  farther  from  the  semiconductor  surface  due  to  the 
insertion  of  phenyl  ring.  Thus,  we  get  a  conclusion  that  the  charge 
recombination  between  injected  electron  and  dye  cation  (/<reci)/ 
redox  couple  in  electrolyte  ( /<reC2 )  of  2  are  both  smaller  than  that  of 
1,  leading  to  a  higher  ^con  of  2  than  that  of  1.  Thus,  combining  with 
their  theoretical  maximum  limit J™ax,  we  believe  that  it  is  the  lower 
77coii  of  dye  1  than  dye  2  that  make  the  Jsc  of  DSSC  based  on  dye  1 
decrease. 

From  the  above  discussion  about  Jsc  and  U0O  it  can  be  clearly 
concluded  that  although  1  has  better  overall  light  harvesting  ability 
(larger  theoretical  maximum  limit J™ax),  the  Jso  \Z0C  and  ri  is  lower 
than  that  of  2  due  to  the  fast  charge  recombination,  in  good 
agreement  with  the  experimental  results  [34].  It  is  the  insertion  of 
phenyl  ring  that  lead  to  the  decrease  of  kYec.  Our  results  also  reveal 
that  the  charge  recombination  process  plays  an  important  role  in 
affecting  the  Jso  Voc  and  77. 

3.3.  Further  optimization  through  molecular  design 

Since  we  have  known  the  main  advantages  and  disadvantages  of 
2  through  the  above  theoretical  analysis  and  the  comparison  with 
the  experimental  results,  these  organic  dyes  (3-5)  shown  in  Fig.  1 
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Fig.  5.  The  optimized  molecular  structures  of  dye-I2  complexes  for  1  and  2. 


were  designed  theoretically,  taking  2  as  a  reference  with  the  aim  to 
further  improve  the  performance  of  the  sensitizer  based  on  the 
following  considerations:  (i)  since  Wang  et  al.  have  demonstrated 
that  the  rigidification  of  bithiophene  is  an  useful  way  to  improve 
the  light  harvesting  ability  by  increasing  the  molecular  conjugation 
degree  [55  , 3  were  designed;  (ii)  through  our  analysis,  we  realized 
that  BTDA  unit  near  the  semiconductor  surface  could  have  a  strong 
interaction  with  I2,  which  could  lead  to  the  fast  charge 

Table  2 

Computed  maximum  absorption  wavelengths  Amax/nrn  (eV)  and  nature  of  the 
transitions  of  1  and  2  corresponding  to  S0  ->  Si  in  dichloromethane  at  TD-CPCM- 
CAM-B3LYP/6-31G(d)  level  with  the  B3LYP/6-31G(d)  optimized  geometries. 


Dye 

State 

Main  configuration 

Amax(e/M  1  cm  !) 

Amax(e/M  1  cm  ^ 

1 

So  — >  Si 

H-l  -*  L  (0.38) 

H  ->  L  (0.52) 

561(2.21)(27,129) 

570(2.18)(18,900) 

2 

Sq  —>  Si 

H-l  ->  L  (0.37) 

H  — ►  L  (0.43) 

501(2.47)(30,178) 

515(2.41)(29,400) 

a  H  represents  HOMO,  L  represents  LUMO,  and  Data  in  parentheses  are  the  main 
configuration  contributions. 

b  The  molar  absorption  coefficients  were  simulated  by  Multiwfn. 
c  Experimental  values  were  obtained  in  dichloromethane  from  Ref.  [34]. 


Fig.  6.  The  absorption  spectra  and  LHE  curves  of  1  and  2  and  their  corresponding J™ax 


LHE  (100%) 


J.-Z.  Zhang  et  al.  /  Journal  of  Power  Sources  267  (2014)  300—308 


307 


Table  3 

Key  parameters  for  deducing  injection  driving  force  (AGinj/eV)  and  the  corre¬ 
sponding  experimental  injection  efficiencies. 


Dye 

Cdye/^V 

Aax 

£dye*AV 

ACinj/eV 

d)  a 

1 

-5.69 

2.21 

-3.48 

0.52 

94.8% 

2 

-5.58 

2.47 

-3.11 

0.89 

94.7% 

a  Experimental  values  were  obtained  in  dichloromethane  from  Ref.  [34]. 


Table  4 

Conduction  band  energy  level  shift  (ACB/eV),  entity  of  charge  transfer  (qCT/e), 
average  electrostatic  potential  on  the  first  layer  of  Ti02  cluster  (VEL/eV),  injection 
driving  force  (AGmj/eV),  theoretical  maximum  limit  of  short  circuit  current  density 
(7™ax/mA  cm  2),  hole-surface  distance  (r/A)  and  position  of  predominant  dye-I2 
interaction  site  for  2—7. 


ACB /eV 

qCT/e 

VE  JeV 

AGmj/eV 

jsr/mA 

cm-2 

r/A 

Position  of 
predominant 
site  (distance/A) 

2 

0.256 

0.36 

-0.21 

0.89 

5.293 

19.9 

2-N2  (6.4) 

2-N1  (14.1) 

3 

0.241 

0.37 

-0.21 

0.86 

6.518 

18.8 

3-N1  (14.2) 

3-N2  (6.7) 

4 

0.257 

0.37 

-0.21 

0.77 

6.000 

18.6 

4-SI  (18.6) 

4-N1  (19.3) 

4-N2  (12.3) 

5 

0.235 

0.37 

-0.21 

0.95 

6.625 

22.1 

5-N2  (6.8) 

5-N3  (7.8) 

5-N1  (14.9) 

6 

0.260 

0.37 

-0.21 

0.81 

8.517 

20.3 

6-N1  (19.4) 

6-N2  (11.0) 

7 

0.240 

0.38 

-0.21 

0.79 

8.918 

20.0 

7-N2  (14.4) 

7-N1  (23.6) 

a  Data  in  parentheses  are  the  corresponding  approximate  distance  between  I2  and 
Ti02  surface. 


recombination  process  between  the  photoinjected  electrons  in  the 
conduction  band  and  the  electrolyte,  thus  4  with  BTDA  unit  far 
away  from  the  surface  were  designed;  (iii)  for  5,  the  molecular 
length  was  increased  by  adding  a  thiophene  ring  aming  to  further 
increase  the  hole-surface  distance  r,  which  is  the  main  reason  that 
contribute  to  the  lower  /<reci  of  2  compared  to  1.  In  a  word,  our  aim 
of  designing  3-5  is  to  further  improve  the  advantages  and  avoid  the 
disadvantages  of  2.  Through  our  calculated  results  shown  in  Table  4 
we  found  the  desired  properties  were  all  presented  in  these  three 
dyes.  However,  none  of  these  three  dyes  are  perfect  due  to  some 
disadvantages.  For  instance,  the  dye  cation  hole-surface  distance  r 
of  dyes  3  and  4  are  18.8  A  and  18.6  A,  which  are  slightly  smaller 
than  that  of  2  (19.9  A),  implying  that  their  /<reci  will  not  be  further 
retarded.  For  dye  5,  although  its  /<reci  has  been  reduced,  there  are 
more  predominant  interaction  sites  than  that  of  2  leading  to  a  faster 
/Creel-  Thus,  for  designing  more  efficient  dyes,  we  integrated  these 
three  molecular  design  principles  into  one  leading  to  the  design  of  6 
and  7  with  just  different  conjugation  order.  The  key  parameters 
affecting  the  performance  of  both  dyes  compared  with  that  of  2 
discussed  above  are  all  collected  in  ^able  4.  From  Table  4,  the  AGinj 
of  dyes  6  and  7  (0.81  eV  and  0.79  eV  respectively)  are  much  larger 
than  0.2  eV,  indicating  that  both  dyes  have  enough  injection  driving 
force  to  inject  electrons  into  the  conduction  band  of  semiconductor. 
Thus,  we  believe  that  both  dyes  have  similar  0jnj  which  are  all  close 
to  one  as  well  as  2.  And  the  corresponding  ACBs  of  6  and  7  are 
0.260  eV  and  0.240  eV,  respectively,  showing  that  dyes  adsorbed 
onto  TiC>2  could  have  noteworthy  effect  on  conduction  band  energy 
level,  leading  to  the  increasing  of  Uoc.  However,  the  difference 
among  ACBs  of  2,  6  and  7  is  so  negligible  that  it  can  not  bring 
significantly  change  on  Voc  in  this  aspect.  The  calculated  J™ax of  6 
and  7  are  32.40  mA  cm-2  and  31.50  mA  cm-2  respectively,  implying 
that  the  theoretical  maximum  limit  J™ax  of  6  and  7  are  much  larger 


than  that  of  2  (22.05  mA  cnrT2).  The  hole-surface  distance  r  sug¬ 
gests  that  the  cation  holes  of  6  and  7  are  slightly  farther  than  that  of 
dye  2  implying  that  the  /<reci  of  both  dyes  are  slightly  smaller.  With 
respect  to  /<rec 2,  although  6  and  7  have  the  same  number  of  pre¬ 
dominant  interaction  sites  with  2,  all  the  predominant  sites  are 
father  than  the  corresponding  ones  in  2.  Thus,  the  /<rec 2  of  6  and  7 
should  also  be  smaller  than  that  of  2.  As  a  consequence  of  the 
above,  it  is  undoubtedly  that  both  dyes  we  designed  will  be  the 
promising  candidates  of  DSSC. 

4.  Conclusions 

The  underlying  origin  of  large  efficiency  difference  for  DSSCs 
based  on  dyes  1  and  2  with  tiny  structure  difference  was  disclosed 
from  the  theoretical  point  of  view.  The  parameters  determining  the 
efficiency  of  the  cell,  including  ACB,  nc,  J™ax,  <t>mj  and  ^con  were 
calculated  through  DFT  and  TDDFT  calculations.  The  results  show 
that  although  2  has  a  little  smaller  theoretical  maximum  limit  J™ax 
than  that  of  1,  which  is  not  beneficial  for  the  enhancement  of  Jso  the 
hole-surface  distance  r  of  2  is  larger  and  the  predominant  inter¬ 
action  site  with  I2  is  farther  from  the  semiconductor  surface  due  to 
the  insertion  of  phenyl  ring,  which  could  lead  to  the  slower  charge 
recombination.  Moreover,  It  gave  us  the  hint  that  larger  hole- 
surface  distance  r  and  the  farther  position  of  BTDA  unit  which 
has  strong  interation  with  the  I2  could  slow  down  the  /<rec  for  the 
DSSC  device.  In  this  context,  two  dyes  derived  from  2  (6  and  7)  with 
better  conjugation  degree,  farther  position  of  BTDA  units  and 
longer  molecular  length  simultaneously  were  designed  to  further 
improve  the  performance  of  DSSCs.  Fortunately,  we  find  that  we 
could  get  dyes  with  our  desired  properties  through  the  reasonable 
rules,  and  6  and  7  are  the  promising  candidates  in  the  DSSC  field 
compared  with  2.  We  hope  our  work  could  provide  a  theoretical 
guidance  for  the  future  research  of  dye  sensitized  solar  cell. 
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